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A purification scheme which utilizes the action of repeated measurements on a (part of a total) 
quantum system is briefly reviewed and is applied to a few simple systems to show how it enables 
us to extract an entangled state as a target pure state. The scheme is rather simple (e.g., we need 
not prepare a specific initial state) and is shown to have wide applicability and flexibility, and is 
able to accomplish both the maximal fidelity and non- vanishing yield. 
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I. INTRODUCTION 

It is well known that in quantum mechanics, the ac- 
tion of measurement affects the dynamics of the system 
just measured in an essential way. This has to be con- 
trasted with the situation in classical mechanics, where 
the effect of measurement can be made as small as one 
wishes. Typical phenomenon reflecting such a peculiarity 
in quantum mechanics has been known under the name 
of "Quantum Zeno Effect (QZE)" Q and has been ex- 
tensively studied recently. It states that if a system is 
frequently measured to confirm that it is in its initial 
state, the change of the state based on the Hamiltonian 
of the system is decelerated, or to state differently, the 
decay of an unstable state is hindered by frequent mea- 
surements. It is widely recognized, however, that there 
is no paradoxical point in such phenomena and the effect 
is solely understood quantum mechanically 0. Indeed, 
the projective measurement is not essential and is just re- 
placed with the generalized spectral decompositions 0, 
and the effect is due to the peculiar short-time dynamics 
of quantum systems, known as the "flat derivative" of the 
survival probability, P(0) — 0. The effect has also been 
examined experimentally and the first announcement of 
its confirmation appeared in an oscillating system Q and 
then in a truly unstable system |^. It is worth while 
mentioning that Raizen's group has observed the devia- 
tion from the exponential law at short times Q and even 
the acceleration of decay when the measurements are not 
frequent enough j^. The latter effect is called the "In- 
verse (or Anti) Zeno Effect (IZE)," which has also been 
studied and discussed ]}. The QZE (and/or IZE) has 
still been explored extensively in the hope of realizing a 
protection scheme for system's coherence against possible 
decoherence by this (or its related) mechanism Q. 

Here another interesting consequence of the peculiarity 
of quantum measurement shall be disclosed. We consider 
a total system that is composed of two (sub) systems A 
and B and measurements shall be repeatedly performed 
only on system A at regular intervals. Our interest lies 



in the (asymptotic) dynamics of system B, which is in- 
directly affected by the measurements on A through its 
interaction with A. It is shown that such indirect mea- 
surements can drive system B to a pure state {purifica- 
tion), irrespectively of its initial state that is mixed in 
general . We can say that the effect of measurement is 
far reaching and profound: It can control even the other 
parts of the system which are not touched directly. In 
the following sections, the mechanism of this purification 
is briefly reviewed (Sec. and is applied to a simple 
qubit system (Sec. IIII|) to show how it works and how it 
can be optimized. Since the entangled state, which is one 
of the key elements in quantum technologies like quan- 
tum computation, information, teleportation, etc. [10|, is 
a pure state of a compound system, this method is used 
to extract entangled states [ill [131 in Sees. II VI and IVl A 
brief summary is presented in Sec. I VII 



II. GENERAL FRAMEWORK OF 
PURIFICATION THROUGH REPEATED 
MEASUREMENTS 

Let the total system consist of two parts, system A 
and system B, and the dynamics be described by the 
total Hamiltonian 



H = Ha + Hi 



Hi, 



(1) 



where Hint stands for the interaction between the two 
(sub)systems. We initially prepare the system in a prod- 
uct state 



Po = 10) (01 ® PS (0) 



(2) 
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at <: = 0. Such a state can be realized, say, if the system 
A is found in the state |0) after the zeroth measurement. 
Notice that system B can be in an arbitrary mixed state 
Pb(0). We perform measurements on A at regular in- 
tervals T to confirm that it is still in the state |0), even 
though the total system A-|-B evolves unitarily in terms 
of the total Hamiltonian H. Since the measurement is 
performed only on system A, the action of such a (pro- 
jective, for simplicity) measurement can be conveniently 
described by the following projection operator 
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o=\^){4>\(g>ii 



(3) 
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Thus the state of system A is set back to |(/)) every after 
T, while that of B just evolves dynamically on the basis 
of the total Hamiltonian H . We repeat the same mea- 
surement, represented by TV times and collect only 
those events in which system A has been found in state 
\(t>) consecutively N times; other events are discarded. 
The state of system B is then described by the density 
matrix 



P^b\n) = (y,(r)) V(0)(vJ(r))~/pM(iV), 
where 

V^{T)^{<j^\e-'"^\cj,) 
is an operator acting on B and 



(4) 



(5) 



Trp 



{V^{r))pBmvl{T)) 



N 



(6) 



is the success probability for these events to occur (yield). 
This normalization factor in 10} reflects the fact that only 
right outcomes are collected in this process. 

In order to examine the asymptotic state of system 
B, consider the spectral decomposition of the operator 
V^{t), which is not hermitian, V0(r) ^ VI(t). We there- 
fore need to set up both the right- and left-eigenvalue 
problems 

V0(t)|u„) = A„|-u„), (w„|V0(t) = A„(w„|. (7) 

The eigenvalue A„ is complex valued in general, but its 
absolute value is bounded [13 



0< |A„| < 1, 



(8) 



which is a reflection of the unitarity of the time evolution 
operator e"*^"^. These eigenvectors are assumed to form 
a complete orthonormal set in the following sense 



\Un){Vn\ = is, {Vn\Um) = 4 



(9) 



Then the operator V^{t) itself is expanded in terms of 
these eigenvectors 



(10) 



It is now easy to see that the A^th power of this operator 
is expressed as 



(11) 



and therefore it is dominated by a single term for large 
A^ 



(12) 



when the largest (in magnitude) eigenvalue Aq is discrete, 
nondegenerate and unique. If these conditions are satis- 
fied, the density operator of system B is driven to a pure 

state 



(13) 



with the probability 



p(^)(Ar) 



large N 



Xor{uo\uo){vo\pBmvo)- (14) 



The pure state |mo), which is nothing but the right- 
eigenvector of the operator (t) belonging to the largest 
(in magnitude) eigenvalue Aq, is thus distilled in system 
B. This is the purification scheme proposed in J^. 

A few comments are in order. First, the final pure 
state |uo) toward which system B is to be driven is de- 
pendent on the choice of the state 10) on which system A 
is projected every after measurement, the measurement 
interval r and the Hamiltonian H , but docs not depend 
on the initial state of system B at all. In this sense, the 
purification is accomplished irrespectively of the initial 
(mixed) state pb(0). Second, as is clear in the above 
exposition, what is crucial in this purification scheme is 
the repetition of one and the same measurement (more 
appropriately, spectral decomposition) and the measure- 
ment interval r need not be very small. It remains to be 
an adjustable parameter. Third, if we can make other 
eigenvalues than Ao much smaller in magnitude 



|A„/Ao| < 1 for n ^ 0, 



(15) 



by adjusting parameters, we will need fewer steps (i.e., 
smaller N) to purify system B. 

It is now evident that the purification can be made 
optimal, if the conditions H15|l and 



|Ao| = l 



(16) 



are satisfied. This condition H16|l assures that we can re- 
peat as many measurements as we wish without running 
the risk of losing the yield (success probability) P'^'^^(A) 
in order to make the fidelity to the target state |uo), 

F<^^\N) = Tib [p^s\n)\uo){uo\/{uo\uo)] , (17) 

higher. Actually, the yield P^"^^ (A^) decays like 



IXol'^iuoWoKvolPBmvo) (18) 



and the condition IjKil) can bring us with the non- 
vanishing yield (uo\uo){vq\pb{0)\vo) even in the N ^ oo 
limit. Therefore the condition Hl()|l makes the two (some- 
times not compatible) demands, i.e., higher fidelity and 
non-vanishing yield, achievable, with fewer steps when 
the condition (|15|) is met. In this sense, the purification 
is considered to be optimal. 
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It would be desirable if an optimal purification can be 
realized by an appropriate choice of the state and/or 
tuning of the measurement interval r and parameters in 
a given system. In the following sections, a few simple 
systems are examined to show how such optimal purifi- 
cations are made possible. 



III. PURIFICATION OF A QUBIT 

As a simplest example, let us consider a total system of 
interacting two qubits. Systems A and B are represented 
by two qubits A and B, respectively, and their two levels 
can be described as the two degrees of freedom of spin- 
1/2 particle. We measure qubit A at regular intervals r 
and examine the state of qubit B, which is in interaction 
with A. The measurement is conveniently parameterized 
as that of a spin- 1/2 particle along a particular direction 
n. After qubit A has been confirmed that its "spin" is up 
along n, its state is projected to the eigenstate of the spin 
operator along n, i.e., \4>) = n-rj^), with t= (ti,T2,T3) 
being the Pauli matrices acting on qubit A. The state 
10) is parameterized in terms of the two angles and Lp 
as 

10) = cos(0/2)e-^'^/^|t) + sin(0/2)e"^/2|i). (19) 
It is an elementary task to write down the projected op- 



erator V^{t) 



-iHt] 



in the following form 



Co 



(20) 



where cr = {ai,<j2, cs) are Pauli matrices acting on qubit 
B and parameters Cg and c = (ci, 02,03) are complex 
valued. Its eigenvalues A± and the corresponding right- 
and left-eigenvectors are easily found 



A± = Co ± c. 



\/c2, c± = ci±ic2, (21) 



{v+\ 
{v-\ 



V2c(c- 


C3) 


1 




V2c(c - 


C3) 


1 




^2c{c - 


C3) 


1 




v/2c(c - 


C3) 



c_|T) + (c-C3)U) 
"(c3-c)|t) + c+U) 

C+(T| + (C-C3)(i| 
(C3-C)(T|+C_(i| 



(22) 



(23) 



Now let the total Hamiltonian of this system be given 

by 

i^=^(l + r3) + ^f(l + .3) 

-I- g{T+(T- + h.c.) + h{T+(T+ + h.c), (24) 

where real parameters g and h are responsible for the 
interaction between the two qubits, A and B. In this 



case, the parameters cq, . . . , C3 in lf^ - (P5|l are expHcitly 
calculated to be 



1 / rOh r9„ 



-(^-sm— + -sm^jcose, (25) 

h . rOh g . Teg\ 
ci = —I \ — sm — - — h — sm — ^ sm cos ip, (26) 

■ I h . T0fi g , t6„ , . „ . 
C2=«(— sm^ — sm^ ) sm0sm(p, (27) 

1 / rOh Te„ 



C3 = 2 I ^ cos ^ ) cos f 



— sm sm — - , (28) 

2\6h 2 6g 2 ' ' ^ ' 



where we have introduced 



LU±=LOA± Ws, 





.UJ+ 




= cos 

2 




sm 

2 


TOg 




. ^Sg 


= COS 




sm — - 


2 




2 



^+4/i2, eg = ^ujl+Ag-^. 

(29) 

In order to illustrate how an optimal purification can 
be achieved in this system, consider the case where we 
measure qubit A along the 3-direction, that is, we choose 
9 = Q and 10) = |t). Then the eigenvalues A± in (|^ and 
the corresponding eigenvectors are 



k+) = IT), (30) 

l"-) = U)- (31) 



Since we have 

sin'ifi, |A_|=.l-^si„=:^, 

(32) 

the purification can be made optimal, e.g., when the pa- 
rameters are adjusted so that hsva.{T9h/2) = is sat- 
isfied. In this case, |A+| = 1 and qubit B is driven 
to a pure state jl), more quickly for larger g satisfying 
sin^(T0g/2) = 1. A similar situation can happen; we can 
extract ||) in qubit B, more quickly for larger h satis- 
fying sin^(T0ft/2) — 1, if we adjust parameters so that 
5sin(T6'g/2) = holds. 

Needless to say, there are cases where such purifica- 
tions are not possible. For example, consider a case where 
we measure qubit A in the 1-2 plane, i.e., — 7r/2. In 
this case, since the parameter co is real, while all the 
other parameters ci, C2 and C3 become pure imaginary, 
the eigenvalues A± = co ± c are degenerated in magnitude 



|A+| = ^/c^ + |cP = |A- 



(33) 



and no purification can occur in this particular case. 
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IV. ENTANGLEMENT DISTILLATION: I 

As is mentioned in the Introduction, since entangle- 
ment is one of the key elements in quantum technolo- 
gies 1 would be useful if the present scheme of pu- 
rification can be used to extract an entangled state as 
a target pure state. Notice that since the target system 
has never been measured directly in the present scheme, 
it is considered to be suited for extraction of such a frag- 
ile pure state as an entangled state. Actually any mea- 
surement on its subsystem that consits of entanglement 
would result in the destruction of the entanglement. 

In order to see an entanglement distillation on the ba- 
sis of the present idea of purification 0, ^| , we consider 
a total system composed of a compound system A-l-B, in 
which an entangled state is to be extracted, and another 
system C. Systems A and B interact with system C sep- 
arately, but do not interact directly with each other. We 
measure system C repeatedly at regular intervals t and 
endeavor to extract an entangled state as a pure state in 
A-t-B. In this section, a simple model, in which all sys- 
tems A, B and C are represented by qubits, is considered 
with a model Hamiltonian 



H 



(1 



(1 



+ g{a^T^ +crfr_ -l-h.c.) 



+ h.c.). 



(34) 



where the Pauli matrices act on system C. It is as- 
sumed here for simplicity that the two systems A and B 
are the same and the Hamiltonian is symmetric under 
the exchange A<->B. 

In order to find the spectral decomposition 1)10(1 of the 
projected operator V^{t) in this case, it turns out to be 
convenient to introduce the Bell states 



1 



= ^ in)±UT) 



|ci>±) = i= |TT)±IU) 



V2 

(35) 

as a complete orthonormal set for A-fB, because we are 
interested in extraction of such entangled states. In- 
deed, the eigenstates of the total Hamiltonian H can 
be found after their classification according to the above 
mentioned A^B symmetry and a "parity" V 
1) A<-!-B symmetric and V = + 



|$+T)\ /n + oj n g + h\ /|$+T)\ 
H\ |$-t) = n n + u; -g + h\ |$-T) , 
l*+i>y \g + h -g + h n J \\^+i)J 

(36) 

2) A<-!-B symmetric and V = ^ 



H 




n n g + h' 

il g ^ h 

g + h g — h + 



/l<f+i>\ 
l<I'-i> , 

Vi*+t); 



3) A<-!-B anti-symmetric and V = — 

H\^-^) = {n + u;)\^-^), 

4) A<-!-B anti-symmetric and V = + 

H\^-l)^L,\^-i). 



(38) 



(39) 



Here l^+t) = |'I'+) <E) |T), etc. Thus the time evolution 
operator e^'^"^ is expressed as 



IT)(TI 



ii)ai 



(40) 



where the summation is taken over the six A<-!-B symmet- 
ric eigenstates of H, denoted as |s), that are given as lin- 
ear combinations of the six states in ((36|l and ((37|l . Owing 
to the A<-!-B symmetry of H, the A<-!-B anti-symmetric 
state Ivt") does not mix with the other (A^B symmet- 
ric) eigenstates. 

We are now in a position to examine the spectrum 



of the projected operator V0(t) = ((/)|e" 
measurement of C projects its state on 



iHT 



a|T)+/3|i), 



If the 



(41) 



the operator reads 



'iEsT 



\mm + \f3\'{i\s){s\i) 



+ a*f3{^\s){s\i) 



(42) 



From this expression, it is evident that the Bell state l^*") 
is always one of the eigenstates of this operator and if the 
measurement interval r is so adjusted that the condition 
LUT = 27r is met, its eigenvalue A,j,- becomes maximum 
in magnitude 



UJT = 2n 



|A*-|-1, 



(43) 



irrespectively of the projected state of system C be- 
cause |ap -I- — 1. This clearly demonstrates a possi- 
bility of entanglement distillation in this simple system, 
in the sense that one of the conditions for optimal pu- 
rification 11()|) can be realized by H43I) and the entangled 
state 1^^) would surely be extracted, only if the other 
eigenvalues of the operator than A^- can be made (much) 
smaller in magnitude. Further details of the analysis of 
such conditions for the entanglement distillation and its 
optimization are found in |lll| and |l2j | in a slightly sim- 



(37) plified case. 
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V. ENTANGLEMENT DISTILLATION: II 

The example in the previous section expUcitly demon- 
strates that we can distill an entangled state in the sys- 
tem A-l-B, through the repeated measurements on the 
other system C that separately interacts with A and B. 
The framework is rather simple and the distillation can 
be made optimal. There is, however, a kind of drawback 
in this scheme. As is clear in its exposition, it is assumed 
that system C, on which the measurement is performed, 
always and simultaneously interacts with both A and B 
and these interactions are crucial for the entanglement 
distillation. Stating differently, systems A and B (and 
C) are not (and/or will not be) able to be separated spa- 
tially, which implies that no entanglement between spa- 
tially separated systems is possible by the scheme pre- 
sented in Sec. lIVI It would not be suited to the situations 
where entanglements among spatially separated systems 
are required, as in quantum teleportation. 

In this section, a resolution to this problem is pre- 
sented. Since the two systems, A and B, an entanglement 
between which is to be driven, are considered to be placed 
at different places, let us consider, instead of system C 
which can no longer interact simultaneously with A and 
B, another quantum system, say X, which is assumed to 
interact with A and B, not simultaneously, but succes- 
sively ^3 • System X plays the role of an "entanglement 
mediator." After such successive interactions with A and 
then B, system X is measured to confirm that it is in a 
certain state. If system X is found in this particular state, 
X is again brought to interaction with A and then with 
B. This process, i.e., X's interaction with A, that with B 
and measurement on X, will be repeated many (N) times 
and we are interested in the asymptotic state of system 
A-l-B in the hope of distilling an entangled state. 

There are a couple of points to be mentioned here. 
First, it is clear that in spite of these modifications, the 
new scheme presented here shares essentially the same 
idea of purification with the previous ones: The dynam- 
ics of the system can be affected, in an essential way, by 
the action of measurement, even if its effect is not direct. 
Second, such a successive interaction would be conve- 
niently treated in terms of a time-dependent (effective) 
Hamiltonian H{t). We may thus avoid possible compli- 
cations caused by the introduction of spatial degrees of 
freedom, still keeping the essential points. 

In order to see how the new scheme works, consider 
again a three-qubit system, A-|-B-f X, for definiteness and 
simplicity. We prepare system X, say in up state |t), 
while the system A-f B can be in an arbitrary mixed state. 
It is assumed that systems A and B are spatially sepa- 
rated and have no contact with each other and that only 
system X can interact with them locally for definite time 
durations. Now consider the following process: 

1. System X is first brought to interaction with system 
A for time duration tj^. The Hamiltonian here is 
given by H{t) = Hq -f- Hxa- Then the interaction 



is switched off and the total system evolves freely 
with the free Hamiltonian i/p for t^. 

2. System X then interacts with system B for time du- 
ration tB^ the dynamics of which is now described 
by another Hamiltonian H{t) = Hq + Hxb- After 
that, the total system again evolves freely with the 
Hamiltonian Hq for tb- 

3. A (projective) measurement is performed on sys- 
tem X to select only up state |t). Other states are 
discarded. 

Then this process is repeated N times: 1-^2-^3— >l—> 

>1^2^3. 

It is shown below that the following choice of the 
Hamiltonians 



i?0 = |(l + a3^) + |(l + CT3^) + |(l + CT3^), 



Hxa = 9A(^i'yt^ Hxb = ffsc^f ct^ 



X„B 



(44) 



actually results in an entanglement distillation in system 
A+B. It is important to notice that the above choice of 
the interaction Hamiltonians is closely connected to the 
details of the process I— >2— >3 and another choice, e.g., 
(Trj^fj^'^'' for HxA(B)i would result not in an entangled 
state, but in a product state, in this particular process. 

The next task is to find the spectral decomposition of 
the projected operator V| defined, in this case, by 



n = (Tie- 



-iHoTB „-i{Ha + HxB)tB 



X e 



iHoTAp-i{Ho+HxA)tA 



IT). (45) 



Since a parity defined hy V = cr^a^ is conserved in this 
system, eigenstates of the operator V\ are easily found. 
Indeed, we can classify every state of system A+B into 
two sectors according to the parity V — and the action 
of the operator V] is closed within each sector. For V — -\- 
states, the action is represented by a matrix A4 



Vf 



ITT)' 
lii) 



_ ^-iuj{tA-\-TA+tB-\-TB) ^ 



ITT)' 
lii) 



(46) 



where its matrix elements read 



Ml2 

M21 



-iLj{tA+2TA + tB+2TB) , 



COS Ca — i sin (a cos 2^^) 



X (cos^B — i sinCs cos2^b), 
sin (a sm2^A sin gBtB, 



-iuj{tB+2TB) 



sing At A sinCs sin 2^b, 



M22 = cosgAtA cosgBtB, (47) 
— states, it is represented by another 



while, for V 
matrix M 



ITi) 
lit) 



_ ^-iuj(tA+tB+'2.TB) 



ITi)' 
liT) 



(48) 
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with its matrix elements 

A/12 = - sin (a sin 2^^ sin Cs sin 2^b , 



^22 



^-iu{tA + 2TA+tB) 



svagAtA sin gstB, 



'iu){tA+2TA) 



COS gAtA (cos Cb — i sin Cs cos 2^^ ) 



(49) 



Here the angles are defined by 

C,A(B) tA{B)^^'^ + 9\b)^ tan2^^(S) = 

(50) 

In order to see the possibility of entanglement distilla- 
tion in this framework, it is enough to consider a much 
more simplified case. Let the two systems A and B be 
treated symmetrically, that is, all parameters are taken 
to be the same for A and B 



9A = 9b — .9) tA =tB = t, TA ^ Tb — T, 
{Ca{B) C; ^A{B) 0- 



(51) 



It is then easy to see that if the parameters satisfy 

cosC — i sin(cos2^ = — e*""^ cosgt, (52) 

an optimal purification of an entangled state \'^) of the 
form 



1*)^^ in) + e^'^UT) 



(53) 



with X = Lu{t + t) is actually possible, provided 

cos gt sin gt 0, uj{t + t) 2mT [n: integer). (54) 

In fact, one can show that |^') is an eigenstate of the 
operator V| 

T/||*) = A^,|*). (55) 
The eigenvalue X^^, is maximum in magnitude 

Avt = -e-3»'^(*+^\ |A*| = 1, (56) 

while all the other eigenvalues remain smaller than unity 
in magnitude, under the conditions (|52|l and H54|) . There- 
fore, we can repeat the process 1^2^3 as many times as 
is required to achieve the desired (high) fidelity, without 
reducing the yield. 

This is an example of (optimal) entanglement distil- 
lations, where the entanglement between two qubit sys- 
tems that are (or can be) spatially separated, is extracted 



through their successive interactions with another qubit, 
on which one and the same measurement is repeated reg- 
ularly. Further details of this model and applications to 
other quantum systems, e.g., extraction of entanglement 
between two cavity modes at a distance, will be reported 
elsewhere. 

VI. SUMMARY 



In this paper, a new purification scheme recently pro- 
posed 9| is applied to a few simple qubit systems to ex- 
plicity show its ability of qubit purification fSec. lIIl|l and 
entanglement distillations fSecs. HVl and rvl) for two-qubit 
systems. The important and essential idea, on which 
these particular examples are based, is to utilize the ef- 
fect caused by the action of measurement on quantum 
systems. 

It should be stressed again that since the basic idea 
is so simple, that is, one has only to repeat one and the 
same measurement without being concerned about the 
preparation of a specific initial (pure) state, this purifi- 
cation scheme is considered to have wide applicability 
and flexibility. Furthermore, it enables us to make the 
two demands — the maximal fidelity and non-vanishing 
yield — compatible. The examples presented in this paper 
just show these characteristics and many variants can be 
devised according to the actual setups. 
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